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Abstract 

Starting from the spectrum of Schrodinger operators on I™, we pro- 
pose a method to detect critical points of the potential. We argue 
semi-classically on the basis of a mathematically rigorous version of 
Gutzwiller's trace formula which expresses spectral statistics in term 
of classical orbits. A critical point of the potential with zero momen- 
tum is an equilibrium of the flow and generates certain singularities 
in the spectrum. Via sharp spectral estimates, this fluctuation in- 
dicates the presence of a critical point and allows to reconstruct 
partially the local shape of the potential. Some generalizations of 
this approach are also proposed. 

keywords : Semi-classical analysis; Schrodinger operators; Equilibri- 
ums in classical mechanics. 



1 Introduction. 

Background. 

Let Ph = —h 2 A + V be a Schrodinger operator where the potential 
V is smooth on E™ and bounded from below. By a standard result in 
spectral theory, P% has a unique self-adjoint realization on a dense subset 
of L 2 (M. n ). As usually, to this quantum operator P^ we can associate a 
classical counterpart with the Hamiltonian function p(x, £) = £ 2 + V(x) on 
the phase space M" x W l . In what follows, we note 4> t the Hamiltonian 
flow of H p = d^p.d x — d x p.d^. 

In the present contribution we are particularly interested in a relation 
between the asymptotic properties of eigenvalues \j(h) of P% : 

P h ^j(x,h) = \j(H)y>j(x,K), 7i— » 0, 
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and the closed orbits of $ t . In geometry spectrum and periodic orbits can 
be related, in a very explicit way, by means of the Selberg and Duistermaat- 
Guillemin 8 trace formulae. In quantum mechanics, the existence of such a 
relation is strongly suggested by the correspondence principle which asserts 
that, in the semiclassical regime h — > 0, certain properties of P^ can be 
related to <&(. In the physic literature, a more precise formulation of this 
principle appeared in the works of Balian&Bloch [T| and Gutzwiller |10j . 
The Gutzwiller formula is usually stated for the trace of the resolvent : 

voi(^) i y ^ 

jaN y ' 7 es E 

where in the r.h.s the sum concerns the closed orbits 7 inside the surface 

S B = {{x, G R n x R n / + V(x) = E}. 

Also VoI(Sb) is the Riemannian volume of He, Sj and A 1 resp. the clas- 
sical action and the stability factor, including the Maslov phase, of 7. 

In mathematics and in physics, such a relation between spectrum and 
periodic orbits provides a powerful tool of analysis and computation. Sec 
e.g. |14) concerning the asymptotic behavior of eigenvectors *&j(x, H) and 
jllj for various applications in quantum chaos. But, for a Schrodinger 
operator on R™, two different divergences occur in Eq.QJ : 

1) The sum over the spectrum can be divergent. If the sum is convergent 
it can also have a divergent behavior when h — ► 0. 

2) The sum over closed orbits is generally divergent. This is the case if |Ay| 
does not decrease fast enough or if the number of periodic orbits of period 
smaller than T is exponentially growing with T. 

Mathematical approach of the Gutzwiller formula. 

As seen above, the question to remove divergences is important and we 
explain below how to proceed. Assume that the spectrum of Ph is discrete 
in some interval [E — e, E + e], a more global sufficient condition for this 
is given in section 2. An accessible problem is to study the asymptotic 
behavior of the spectral distributions : 

1 {E,h,y)= <P( Xj{h) h ~ E ), asft->0, (2) 

|A 3 (ft)-£|<£ 

where tp is a function chosen to remove the divergences. To justify the 
terminology, observe that the truncated spectral distribution : 

\\j{h)-E\<E 
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acting on a function ip shifted by E and scaled w.r.t. h provides : 



In general, it is not possible to compute the spectrum of P% and one 
motivation is to derive statistics about eigenvalues. For example, in Eq.J2J 
the formal choice of ip as the characteristic function of [—77,77], < 77 < e, 
determines the number N(h) of bound states in [E — rjh, E + r/h]. Under 
certain conditions, it can be proven that N(h) is proportional to h 1 ^ n 
(Weyl-law) . Accordingly, for n > 1 this implies that the finite sum defining 
"f(E, h, tp) involves a large number of eigenvalues as h — > 0. 

A second aspect is that the asymptotic expansion of j(E, h, ip) involves 
the classical dynamics in a very explicit way. We recall that E is regular 
if Vp(x,£ l ) on and critical otherwise, a critical point (xq,£o) of 
p is a fixed point of $ t since H p (xo,£,o) — 0. When E is not critical the 
asymptotics of Eq.J2J is well determined by the closed orbits of $t on E^. 
For the full treatment of this problem, and a complete formulation of the 
asymptotic expansion, we refer to . 

We explain now why the problem stated in Eq.JSJ) leads to a mathe- 
matically rigorous version of the Gutzwiller formula. First, for each h > 
the sum is finite and a fortiori convergent. A convenient choice of ip also 
ensures that this quantity has an asymptotic expansion when h — > inde- 
pendently from the choice of e. On the other side, it will be proven that 
only the periods of $ t inside supp(^3), the support of the Fourier transform 
of ip : 



contribute in the asymptotic expansion. This principle is useful since when 
supp((y3) is compact then finitely many closed orbits of contribute and 
the second divergence is solved. Hence if <p € (R) , the space of smooth 
functions with compact support, ip is in the Schwartz space 5(R). Since ele- 
ments of iS(R) are smooth with exponential decay at infinity, no divergence 
occurs and the size of e is irrelevant in the semi-classical approximation. 

Finally, in Eq. j2J) the scaling w.r.t. h is important. With this choice and 
via Fourier transform considerations, we can use the propagator Unit) = 
expiitPft/h), solution of the Schrodinger equation : 



to obtain a precise control w.r.t. h. Roughly, U%(t) can be expanded w.r.t. 
h via a so-called WKB approximation. This expansion also provides the 
explicit relation with the classical dynamics. The precise technical justifi- 
cations are given in section 3. 





(3) 



ihd t U h (t) = P h U n (t), 
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Critical values and contributions of equilibrium. 

We have outlined the heuristic relation : 



lim 7 (E, ft, if) ^ {(t, i,()e!xS £ / $ t (x, = (a;, £)}• 



In the r.h.s any point of a periodic orbit appears only at times kT, 

k G Z, where T is the primitive period of the orbit. But an equilibrium 
point (xo,£o) satisfies $4(20, £o) = ( x o,£,o) f° r all t- Hence when E is no 
more a regular value the nature of the set of fixed point changes and some 
new contributions appear in the asymptotic expansion. 

When E = E c is critical, the asymptotic behavior of Eq.J2J) is more 
complicated and is closely related to the geometry of the flow inside ^e c - 
For a non-degenerate critical point, i.e. d 2 p is an invertible matrix when 
dp = 0, the reader can consult The problem is treated there for quite 
general operators, also including the case of a manifold of critical points, 
but for supp(</3) small around the origin. For Schrodinger operators on R n 
and supp(<y3) compact but arbitrary, the results of |2] are improved in |13| . 

Two important problems occur in presence of critical points. First, at 
every point where dp = the surface E£ c and the metric of S^ c are not 
smooth. Next, the determination of the asymptotic expansion w.r.t. ft can 
be very difficult. The point is that j(E, ft, <p) can be expressed in terms of 
oscillatory integrals : 



Note this oscillating factor ft -1 , precisely imposed by the scaling w.r.t. ft 
in Eq. (J2J). Via the WKB approximation, the phase / is related to the flow 
so that the asymptotic behavior of 1(H) is determined by the closed orbits. 
The technical problem is that, in presence of an equilibrium, / has some 
degenerate critical points. The stationary phase method cannot be applied 
and the asymptotic expansion of 1(h) is radically different : e.g. some 
terms ft Q , a G Q and powers of log(ft) generally appear in this setting. 

Results and strategy. 

Our objective is to relate some variations in the discrete spectrum of P% 
with the presence of fixed points for the classical system. Conversely, an 
attempt is made to prove that the knowledge of such a spectral fluctuation 
can describe the singularity of the potential. In theory, such a determina- 
tion is possible since the contributions of equilibriums are highly sensitive 
to the local shape of V . 

We will consider the case of a potential V with finitely many critical 
points x J attached to local homogeneous extremum of V. An immediate 



/(ft) 




RxK 2 ' 
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consequence is that p admits, locally, a unique critical point (xq,0) on the 
surface £ £J = {(x,£) € K 2n / £ 2 + F(x) = V(a4)}< A typical example is a 
polynomial double well in dimension 1 where 3 critical points occur at the 
2 minima and at the maximum of V. 

In fact, starting from a more precise relation : 

lim 7 (.E, h, ip) ^ {(*, x, f) G supp(<£) x S B / $ t (x, = (a, 0}, 
the core of the proof lies in 2 facts : 

- Equilibriums have a continuous contribution w.r.t. the time t. 

- A convenient choice of supp(^) erases all other contributions. 

Here, 'continuous contribution w.r.t. V means that a fixed point contribute 
to the asymptotic expansion of -f(E J c ,h,ip) in the form h a log(ft)' 3 (D, ip) 
with supp(£>) = R. Contrary to standard periodic orbits whom contribu- 
tions are supported in the set of periods, such a term supported on the 
line cannot be erased just by shrinking the support of ip. For example, if 
supp(<£) contains no period of the flow the analysis easily follows if we view 
j(E, H, p) as a function of E : 

- The order w.r.t h of j(E, h, ip) changes when E — > E 3 C . This indicates the 
presence of an equilibrium for $ t , a fortiori of a critical point for V. 

- This discontinuity at E 3 C describes the shape of V. 

2 Hypotheses and main result. 

Let p(x, £) = £ 2 + V(x), where the potential V is real valued and smooth 
on R™. To this Hamiltonian is attached the ^-differential operator Pn = 
—h 2 A + V(x) and by a classical result Pn is essentially self-adjoint if V is 
bounded from below. 

Remark 1 We are here mainly interested in the case of Schrodinger op- 
erators but a generalization to an ft-admissible operator (e.g. in the sense 
of |17|) of principal symbol £ 2 + V(x) is given in the last section. 

First, to obtain a well defined spectral problem, we use : 

(Hi) V G C 00 (M"). There exists C £ R such that liminf V > C. 

oo 

Note that (Hi) is always satisfied if V goes to infinity at infinity. Now, 
consider an energy interval J = [Ei,E2\ with E^ < liminf!^. In the 

oo 

following we note : 

J(e) = [Ei-e,E 2 +e}. (4) 
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For e < £o the set p _1 (J(e)) is compact. By Theorem 3.13 of ^7] the 
spectrum a(Pn)r\J(e) is discrete and consists in a sequence Ai(fi) < Aa(?i) < 
... < Xj(h) of eigenvalues of finite multiplicities, if e and h are small enough. 
The main tool of this work will be the spectral distribution : 

y(E,H,<p)= £ (5) 

more precisely, the asymptotic information contained in this object. To 
avoid any problem of convergence we impose the condition : 

(H2) We have tp £ Cq°(R) with a sufficiently small support near the origin. 

Remark 2 (H2) is used to erase contributions of non-trivial closed orbits 
and can be relaxed to <p 6 C^°(R) with a weaker result. A more precise 
description of supp(<£>) is given in Lemma jlOl For a non-degenerate min- 
imum, it is more comfortable to assume that supp(<£>) contains no period 
of d$t(z ). Some singularities, not relevant here, are generated by these 
periods and we refer to ^31 for a detailed study of these contributions. 

To simplify notations we write z = (x,£) € R 2 " and = p _1 ({i?}) and 
we use the subscript E c to distinguish out critical values of p. Of course 
one can also work with T*R™ ~ R" x R". In J there is finitely many critical 
values Ec,...,E l c and in p _1 ( J) finitely many fixed points Zq, z™, m > I. 
We impose now the type of singularity : 

(H3) On each H E j the symbol p has isolated critical points Zq = (cCq,0). 
These critical points can be degenerate but are associated to a local ex- 
tremum of V : 

V(x) =E C + V 2k (x) + 0(\\x - x 3 Q \\ 2k+1 ), k e N*, (6) 
where Vik, homogeneous of degree 2k, is definite positive or negative. 

Remark 3 For non-degenerate singularities we can apply the results of |^j 
and the extremum condition is not really necessary. 

The next assumption, erases the mean values in the trace formula : 
(H 4 ) ip is flat at 0, i.e. ipW{0) = 0, Vj £ N. 

We could weaken (H4) to (p^(0) — 0, Vj < jo, where jo depends only 
on the degree of the singularities of V (cf. section 4). Note that such a ip 
exists. Pick g e Cg°(R), supp(g) C [-M, M], then ip(t) = t 2 ^g(t) satisfies 
our hypotheses. In this case, we can pick g even so that ip is real. 
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Finally, to relax a bit (H2) we need a control on the contribution of 
closed orbits. To do so, we impose the classical condition : 

(H5) All periodic trajectories of the flow are non-degenerate. 

Non-degenerate closed orbits are those whose Poincare map does not admit 
1 as eigenvalue and are isolated. The main result is : 

Theorem 4 Assume (-Hi) to (H4) satisfied. As H tends to + , we have : 

, h (<D(h°°) if se [E 1 ,E 2 ]\{El,...,E l c }, 
7ls ' ,ip) ~ I o{f 3 {h)) ils = Ei je{i,...,i}, 

where each fj(h) has a finite order w.r.t. h. 

Precisely, if s = E 3 C carries a single minimum of degree 2k we obtain : 

fj(K) = C(n,fc s (7) 

But for a local maximum of V we can obtain a logarithm of h : 

fj(H) = C(n,k,^+^- n \og(hy, j = or 1. (8) 

In fact if the critical surface carries more than one critical point then fj is 
the sum of their respective contributions. Note that for n = 1 and k > 1 
the singular term has negative order w.r.t. h. A more detailed formulation 
of each fj(H) is given in Propositions ll2ll31 An interesting property is that 
the singularity of j(s, h, ip) describes partially the singularity of V. 

Corollary 5 Assume that T*e c carries exactly one singularity (xq,0). Then 
the discontinuity of 7(3, h, ip) at s = E c determines the degree of the critical 
point and the spherical average of the germ of V in xq. 

This principle is limited in presence of multiple equilibriums on the same 
surface since the sum of contributions of each critical point could lead to a 
compensation. In (-H2) the condition that supp((^) is small implies a very 
accurate spectral estimate (e.g. by a Paley- Wiener estimates for ip). It is 
possible to relax this assumption but the result is weaker : 

Corollary 6 Assume (Hi), (-H3), (^4) and (H5) satisfied and that tp £ 
Co°(M), then we obtain : 

7 (s, h, ^) = 0(1) if s e [E ll E 2 ]\{E 1 c , ...,E l c }. 

For critical values of p, estimates are the same as in Theorem 0] 
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The justification, given in section 4, is that in this case the asymptotics 
is given by a finite sum over periodic orbits of energy s. This result is 
weak if the singularity of V is non-degenerate since the equilibrium has a 
contribution of degree w.r.t. h. (cf. Propositions ll2ll31 or section 3 of [2]). 
Finally, we would like to emphasize that a maximum is more difficult to 
detect contrary to a local minimum which is isolated on the energy surface. 



3 Oscillatory representation. 

The construction below is more or less classical and will be sketchy. The 
only change is the choice of a more global localization around J = [Ex, E 2 \- 
Strictly speaking, with (Hi), we could also consider ] — oo,E 2 ]- Let be 
ip E S(M) with if e Cg°(R), we recall that : 

1 {E,h, V )= Y, y( Aj( ^~ E ), J{e) = [Ex -s,E 2 + e ], 

X:j(H)eJ(e) 

with p~ 1 (J(e)) compact in T*W l . For e > small enough, we localize 
around J with a cut-off O £ Cq°(]Ei — e,E 2 + e[), such that 9 = 1 on J 
and < 9 < 1 on R. We accordingly split-up our spectral distribution as : 



with 



j(E,h,<p) = 7i (£,/(,,</?) + ~f 2 (E, h,<p), 

A 3 (fe)eJ(e) 

Aj-(ft)eJ(E) " 
Since cp 6 S(M.) a classical estimate, see e.g. Lemma 1 of 0], is : 

-yi(E, H, <p) = 0(h°°), as ft ^ 0+ (9) 
By inversion of the Fourier transform we have : 

R 

The trace of the left hand-side is 72 (E, K, cp) and Eq. @ provides : 

*y(E,h,ip) = ^Tr / e^0(t)exp(-|p fi )e(P ft )rft + 0(tf°). (10) 

Z7T J n 
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Eq. I|1U[) is very close to the classical Poisson summation formula on S 1 since 
the r.h.s. is expressed below in term of the classical dynamics and this 
asymptotic relation justifies the terminology of 'trace formula'. Moreover, 
this formulation shows that the scaling w.r.t. h, imposed in the definition of 
"f(E, h, ip), is the best one since we will solve the semi-classical propagator 
homogeneously w.r.t. h. 

Let Unit) = exp(— jj-Ph) be the quantum propagator. We approximate 
Un(t)Q(Pn) by a Fourier integral operator (FIO) depending on %. Let A be 
the Lagrangian manifold associated to the flow of p : 

A = {(t,T,x,{,y,ri) e T*R x T*R" x T*R n : T=p(x,t), (x,0 = 3> t (y,r,)}, 

and /(R 2n+1 , A) the class of oscillatory integrals based on R 2 ™ +1 and whose 
Lagrangian manifold is A. The next result is a semi-classical version of a 
well known result on the propagator, see e.g. Duistermaat 7 . 

Theorem 7 The operator Un(i)Q(Ph) is an H-FIO associated to A. For 
each JVeN there exists U^^it) with integral kernel in Hormander's class 

7(R 2rl+1 , A) and R^\t) bounded, with a L 2 -norm uniformly bounded for 
< h < 1 and t in a compact subset of R, such that : 

U h (t)Q(P h ) = ui N i(t) + h N R[ N) (t). 

This result provides the existence of an asymptotic expansion in power of h 
with a remainder that can be controlled since supp(<p) is a compact. After 
perhaps a reduction of e, this remainder (t) is estimated via : 

Corollary 8 Let Q 1 e Qj°(R), with 8 X = 1 on supp(6) and supp(9i) C 
}E 1 - 2e, E 2 + 2e[, then VA^ £ N : 

Tr(6(P ft M^^)) = i-Tr J 0(t)ei tE U^(t)e 1 (P h )dt + 0{h N - n ). 

For a proof of this result, based on the cyclicity of the trace and a priori 
estimates on the spectral projectors (see we refer to 0|. For the 

particular case of a Schrodinger operator the BKW ansatz shows that the 
integral kernel of Uq} it) can be recursively constructed as : 

K { h N \t,x,y) = — IL; / bi N \t,x,y,OeU s ^-^M, 

b]p =b + m 1 + ... + h N b N , 
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where S satisfies the Hamilton- Jacobi equation : 

p(x,d x s(t,x,t)) + atS(t,x,o = o, 

with initial condition S(0, x,£) — ( x t0- I n particular we obtain that : 

{(t,d t S(t,x,ri),x,d x S(t,x,ri),d ri S(t,x,ri), -77)} C A, 

and that the function S is a generating function of the flow, i.e. : 

Md v S(t, x, r,), 77) = (x, 8 x S(t, x, r?)). (11) 

We insert this approximation in Ea. l|10|l . we set x = y and we integrate 
w.r.t. x. Modulo an error (D(h N ~ n ), we obtain that j(E, h, ip) equals : 



1 



(2nh) n 

IxT 



>M s ^-^+ tE )ai N \t,x,Omdtdx(%, (12) 



where a^\t, x, rf) = b^(t, x, x, rf). 

Remark 9 By Theorem 3.11 & Remark 3.14 of J7|, Q{Pn) is h- admissible. 
Moreover, the symbol is compactly supported in p~ 1 ([Ei — e, E 2 + s]). This 
point allows to consider only oscillatory integrals with compact support for 
the evaluation of the spectral distributions. □ 

Microlocalization of the trace. 

If ip G C^(T*R™), we recall that ip w (x, HD X ) is the linear operator obtained 
by Weyl-quantization of ip, he. : 

^( x ,HD x )f(x) = -±— J e*(*-y<Z^(^±y,£)f(y)dyd(. 

Mainly, the contribution of an equilibrium zq G Sb c can be reached via : 
lzo (E c ,h,<p) = ^Tr I e l ^<p(t)4> w (x,hD x )eM-fPhMPH)dt, (13) 



where ip G C^°(T*K. ra ) is equal to 1 near zq. This principle will also be 
useful to obtain a weak generalization in presence of multiple equilibriums. 

We recall some basics on symbolic calculus with FIO. Hormander's class 
of distributions with Lagrangian manifold A over R n is noted J(R™ , A) . If 
(xo,£o) G A and ip(x,9) G C°°(M. n x R^) parameterizes A in a sufficiently 
small neighborhood U of (a;o,Co), then for each u% G 7(E™,A) and x G 
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Co°(T*R n ), supp(x) c U, there exists a sequence of amplitudes Cj(x,9) £ 
Cg°{R n x R N ) such that for all L E N : 

X^(x,^K- ^ h j I( Cj e^)+0(H L ). 

-d<j<L 

Hence, for each N E N* and modulo an error (D(h N ~ d ), the localized trace 
7 Z0 (£' C , fi, yj) of Ea. H13|) can be written as : 



SxR 2 " 



To get the right power — d of h, we apply results of Duistermaat [7] on 
the order : fi,-pseudo-differential operators ip w (x, HD X ) are of order w.r.t. 
1/h. Since the order of Un(t)Q(Pn) is — t, we have : 



Multiplying by ip(t)eh tEa and passing to the trace we find Ea. (|14|l with 
d = n and we write again ai (t, cc, 77) for oi (t, cc, n, 77). In particular : 

4 0) (t,x,x,77) = i/j(x,ri)ao(t,x,x,r)), (15) 
is independent of fi. and is compactly supported w.r.t. (x, rf). 



4 Proof of the main result. 

Classical dynamics near the equilibrium. 

A generic critical points of the phase function of Eq. i|12|) satisfies : 



E=-d t S(t,x,0, 
x = d^S(t, x, £), ^> 

C = d x s(t,x,o, 



p(x,0=E, 



where the right hand side defines a closed trajectory of the flow inside Eg. 
Note that equilibriums are also included. By the non-stationary phase 
lemma, outside of the critical set the contribution is 0(h°°). 

Let be E c any critical value in [E\ , E2} and zq an equilibrium of Sb c . 
We choose a function ip E Cfi° (T*R n ) , with ip = 1 near 20, hence : 

■y 2 (E c ,h,<p) = -^Tr / e^^t^a:, HD x )exp(-UP h )Q{P h )dt 

Z7T J /l 

+ ^Tr / e^^)(l-V tu (^,^))cxp(-^)e(P ft )*. 
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If there is no other singularity on £# c with (H5) the asymptotic expansion 
of the second term is given by the semi-classical trace formula on a regular 
level. For finitely many critical point on £_e c , we can repeat the procedure. 
The first term is micro-local and precisely generate the singularity in Theo- 
rem We note the discrete set of critical points z 3 in J). The next 
result provides a global information on the periods of the classical flow. 

Lemma 10 There exists a T > 0, depending only on V and J = [E\, E 2 ], 
such that $ t (z) ^ z for all z e p-^J^fi and all t e] - T,0[U]0,T[. 

Proof. If flp is our hamiltonian vector field and z — (x, £) we have : 

||flp(*i) - H p (z 2 )\\ 2 = 4||a - 6iP + ||V B V(a:i) - V^(x 2 )|| 2 . 

When z\ and z 2 are in the compact j? _1 (J) there exists b > such that : 

||V x y(x 1 )-V x y(^)|| <6||si-x a ||. 

Hence, there exists a > such that : 

H-ffpfzi) - i? p (^ 2 )| I < a\\zi - z 2 \\, Vzi, z 2 £ P^i-J)- 

The main result of JS] shows that any periodic orbit inside p _1 (J) has a 
period r > 2?r/a > 0. The lemma follows with T := T(V, J) = 2-nja. ■ 

Remark 11 The result of ^Hl is optimal (harmonic oscillator). Note that 
T is decreasing if one increase the size of J. Lemma ^] provides a total 
control on the r.h.s. of the trace formula : if (p G C^°(] — T,T[), the only 
contribution arises from the set {(t,Zo), t £ supp(t^)}. 

Now, we restrict our attention to the singular contribution generated by one 
critical point. As pointed out in section 2, for a non degenerate extremum a 
minor technical problem could occur. We recall that the linearized flow d$ t 
is the differential of the flow $ t w.r.t. initial conditions z = (x,£). When 
zq is a critical point of p, the linear map z <— * d$ t (zo)z can be interpreted 
as the Hamiltonian flow of z 1— > (d 2 p{zo)z 1 z). After perhaps a change of 
local coordinates near xq, we can assume that d 2 V(xa) is diagonal. If xq 
is a maximum of the potential d^t(zo) has no non-zero period which ends 
immediately the discussion. If xo is a minimum d<5> t {zo) is elliptic with 
primitive periods (Ti, ..,T„) generated by the eigenvalues of d 2 V(x ). But 
the constant b of Lemma ITU1 is certainly bigger than the spectral radius of 
cPV(xo) and hence we have the inequality T < min{Ti, .., T n }. Following 
the approach of or if supp(<^) c] — T, T[ the associated contribution 
is smooth on supp(^)\{0}. For a degenerate critical point as in (H3) a 
surprising result, established in 0E1j is that the only singularity is located 
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at t = 0. Hence no extra assumption on ip is required. 
The trace as an energy function. 

As seen in section 2 it suffices to study the localized problem : 



Here ip G C^°(T*R™) is micro-locally supported near z (cf section 2). For 
the convenience of the reader we recall the contributions of equilibriums in 
the trace formula. We note S(§" -1 ) the surface of §™ _1 and in the next two 
propositions it is understood that conditions (Hi) to (H 3 ) are satisfied. 

Proposition 12 If Xo is a local minimum we have : 

j,ieN 2 

where the Ajj are some distributions. The leading coefficient is : 

fi9+ft-" S fl* n ) f \V2k(ri)\-*dri [ v(u 2 + v^u^v^dudv. 
(27rj™ J J 

S™- 1 R+XR+ 

Proposition 13 If xo is a local maximum we have : 

Jzo (E c ,h,<p)~h%+%- n ]T E ^ + ^og{h) m A^ m {y). 
m =°. 1 i,'eN 2 



If =^ i N, the first non-zero coefficient is given by : 

+ (T „, fc ^} -L-^i y |V 2fc (»?)|-*d»7- 

The distributions T n! fc are respectively given by : 

(T n ,kM= j(Ctk\^~ 1 +C-^t\ n _ ! ^- 1 Mt)dt, if n is odd, 

R 

(T n ,k,V>)=C- k J if ^ ^ 



is even. 



But if e N and n is odd then the top-order term is : 

c n , fc iog(ft)ft9+ft- n ^|^ y iv^ fc (*?)i-*dT? y i*r^-v(*)*- 
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Finally, if ^ 6 N and n is even, C^ k = C n k and we have : 

S"- 1 R 

Remark 14 To emphasize the consistency of these results we precise that 
C n ,k, &nk are n on-zero universal constants depending only on n and k. 
Such terms H a and h a log(fi), a s Q never appear if E is regular. 

For a proof we refer to |3J|S] resp. for a local minimum and maximum. The 
case k = 1 was treated in 0. With (H$) and E 1 regular, we have : 

fcl-n 00 

j(E,h,^) ~ — LVol(E £ )^(0) + 5> 1 - B+ ^(0)(O) 

OO 

We refer to ^Hl for a proof. In the r.h.s. the sum concerns periodic orbits 
p of energy E and is finite since supp(<p) is compact. Here S p , p p and T p 
are resp. the action, the Maslov-index and the period of the closed orbit p 
and both Cj, D p j are differential operators of order j. If <p satisfies (H4) 
we have Cj(<p)(0) = and for each s £ [Ei, E 2 ] regular : 

oc 

7 (s,h,<p) ~ E e^V^E^^)^'- (16) 

pes s 3=0 

We accordingly obtain that : 

7(jj, ft,p) = 0(1), VaG^i.^M^ 1 ,...,^}. (17) 

This point will justify Corollary By Lemma 1101 we have T p > T uni- 
formly w.r.t. s 6 [EijEz]. Hence if s is not critical and (#2) is satisfied 
the sum over the periods of Eq. (|16f) is simply and in Eq. (|17l) we obtain 
in fact 0(h°°). Note that (#5) is not required here. For s — E™ critical 
there is a continuous contribution w.r.t. t in the spectral distribution. A 
fortiori, a choice of ip flat at the origin does not erase this term. We have : 

7 (£ c "\^)~E/,m 

3=1 

where N m is the number of equilibrium on and each fj(ti) is given by 
the leading term of Propositions 1121151 ■ 
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Note that the bottom of a symmetric double well gives a similar answer as 
a single well of same nature. Hence without microlocal considerations it is 
difficult to distinguish these 2 different settings. 

Proof of corollary\^ First, the Weyl-law for regular energies : 

j(E,h,<p) ~ (2^) 1 -"< ( 5(0)Lvol(E £ ), 

computes the dimension n. Now assume given a critical value E c with a 
single critical point. The only choice of the spectral function ip allows to 
detect E c via the singularity f(h) of Theorem QJ The knowledge of f(H) 
determines the order of the contribution. For example, if : 

f(h)~ch a io g (h), 

the critical point is a maximum and a computes the degree 2k of the 
singularity. With ip, the knowledge of k allows to compute the quantity : 

R 

A fortiori C determines the average of |V2fe| - ^ on S n_1 . Without \og(h), 
the nature of the critical point can be detected by a symmetry argument 
w.r.t. ip since we a priori know n and k. In view of Propositions 1 1 21 1 31 we 
can choose ip odd, even, symmetric or non-symmetric w.r.t. the origin to 
conclude. Note that if (p is not even ip is a priori complex valued. I 

The spherical average of V^k is a Jacobian. For example we have : 
Je-\ v *>Wdx = ±-r(±) J \V 2k (rj)\-^drj. 

A similar result holds for the pullback f(V2k(%)), if / G r^ _1 dr). 

Remark 15 Enlarging the list of singularities would provide a bigger "dic- 
tionary" . The case of non-homogeneous singularities for V is still an open 
problem, in particular because the determination of an explicit asymptotic 
expansion w.r.t. h can be very difficult. 

We propose now 2 slight generalizations of the main result, 
a) Effect of a sub-principal symbol. 

Because of some recent developments of Helffer&Sjostrand for Witten Lapla- 
cians, see e.g. |12j for an overview and references, we show shortly how to 
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extend the result of Theorem0]to the case of an ^.-admissible operator. For 
example, the Witten Laplacian on zero-forms is : 

Afl = -ft 2 A + \\Vf{x)\ 2 - J A/(:r), / e C°°(R»), 

whose symbol p(x,£) = Pq{x,Q + ftpi(cc,£) depends on ft. More generally, 
it is possible to consider operators P% of symbol p^ ~ ^Pj (Borel sum) 
with principal symbol po(x,£) = £ 2 + V(x) an d a subprincipal symbol 
Pl 7^ 0. Starting from the results of section 3 we proceed as follows. 

To each element u% of J(M n , A) we can associate canonically a principal 
symbol s a pi i nc (uh) , where S is a function on A such that £dx = dS on 
A. In fact, if un can locally be represented by an oscillatory integral with 
amplitude a and phase ip, then we have S — S v — (p o i^ 1 and cr pr i nc (uf l ) 
is a section of \A\i tg> M(A), where M(A) is the Maslov vector-bundle of 
A and |A|a the bundle of half-densities on A. When p\ ^ 0, in the global 
coordinates (t,y,rj) on A, the half-density of Un(t) is given by : 

t 

u(t,y,rf) = exp(i J pi($ s (y, -r]))ds)\dtdydr]\i . (18) 
o 

For this expression, related to the resolution of the first transport equation 
for the propagator, we refer to Duistermaat and Hormander Accord- 
ingly, the FIO approximating the propagator has the amplitude : 

t 

a(t, z) = a(t, z) exp(i J p\{<& s {z))ds) . 
o 

Since zq is an equilibrium we have p\(^ a (zo)) = Pi{zq), Vs, and : 

~a(t,z ) = p(ty tMzo) - (19) 

If the subprincipal symbol vanishes at the critical point, which is the case 
in a lot of practical situations, the trace formula remains the same. If 
Pi{zo) 7^ 0, by Fourier inversion formula we replace ip(t) by (p(t + pi(zq)) 
in all integral formulae of Propositions 1121131 Note that, with (H4), this 
has absolutely no effect for the mean values and hence on the main result. 

b) A micro-local approach. 

We inspect now the case of an energy surface supporting more than one crit- 
ical point, but with a much more restrictive method. Let be K = p~ 1 ( J) C 
T*M. n and ro = \ inf d(zi, Zj), where d is any distance on T*R n . Each open 

ball B(z,tq) C T*R" contains at most 1 critical point for each z G K. 
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Clearly, we can cover a compact neighborhood of K by a finite number 
of balls B(z,ro). With a partition of unity, adapted to this covering, we 
obtain : 

N 

J2 ^(s, TiDx) = W, on C™(K). 

3=1 

For each s € J, we obtain : 

Tr J 0(t)Q{P h )e^- s) dt = ^ Tr / 0{t)^{x, hD x )Q{P h )e^- s Ut. 

E J = 1 R 

Note that the r.h.s. is studied in section 2. By the same argument as 
before, if E s n supp(V>j) contains no critical point we obtain : 

k 

And if there is exactly one critical point zq G E s in swpp(ipj) we have : 
Tr J £(t)i>J(x,hD x )e(P h )eT^ p *- s Ut = ^(z )f 3 (h), 

K 

and by construction no cancellation can occur. 

Remark 16 In CorollaryEl we have considered for the flow. A similar 
result holds for a chaotic dynamics and an isolated degenerate closed orbit 
can be treated as in JSj. Finally, using the results of ^S] one can extend 
Corollary to the case of families of periodic orbits of dimension d < n. 
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